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Abstract 

K* \ Accumulation and annihilation of weakly interacting massive particles (WIMPs) in the earth 

^^ ' and the sun may be observed by the resulting neutrino signal. We demonstrate that, for 
certain parts of the SUSY parameter space, present and future neutrino telescopes can expect 

Qv ' a higher rate of events due to tau-neutrinos than of events due to muon-neutrinos. We 

f^ \ show how Vt <-> j/^ oscillations, and in the case of the sun also scattering, absorption and 

^^ ' regeneration, modify the signal expectation. We find that most currently proposed neutralino 

^p \ models predict equal or increased muon fluxes if Vr <^ ^a oscillations are included. 



X 



' , 1 Introduction 

Mh. 

(D 
,^ . Non-baryonic dark matter is believed to exist in order to explain the available data on large 

scale structures in the universe ||l|. A good candidate for such non-baryonic dark matter is the 

lightest supersymmetric particle (LSP) and in particular the lightest neutralino, xo- Cosmolog- 

ical arguments restrict the mass of the neutralino to be lighter than a few TeV [g] and combined 

5^ \ with data from accelerators set lower limits of about 50 GeV 1^]. WIMPs populating the galactic 

halo would loose energy by scattering in the sun and the earth and if slow enough they would 

be trapped by the gravitational potential. These trapped WIMPs then annihilate with each 

other to produce quarks, leptons, gauge or higgs bosons which then decay and thereby produce 

neutrinos. The differential energy flux of these neutrinos is given by 

dNi r ^ dNF,^ ,^. 

-dE = ^^^\^^^E- ^^) 

t 

where T is the annihilation rate within the sun or earth, R is the distance to the detector location, 
Bp is the branching ratio for WIMP annihilation into the different final states F and ^^'' is 
the differential number of neutrinos of flavor i produced in the decay of F. Detailed calculations 
of observable event rates exist for a wide range of parameters of the Minimal Supersymmetric 
Standard Model (MSSM) [Q, |5|. The evaluated signal typically consists of neutrino-induced 
muons which, due to their long range, have an experimentally promising signature. 

'Permanent address 



Here we want to draw attention to the tau-neutrino for the following reason. Annihilation of 
WIMPs can be a significant source of tau-neutrinos resulting in a flux of tau-neutrinos which is 
up to 6 times higher than the flux of muon- or electron-neutrinos. This is realized for those parts 
of the allowed MSSM parameter space in which the branching ratio for annihilation into a pair of 
taus is large. The tau decay yields always a tau-neutrino and in 17% also an electron- or muon- 
neutrino. A high tau-neutrino flux has strong implications since present data of atmospheric 
neutrinos provide evidence for strong Vj- •s-^ v^ mixing |p. Tau-neutrinos would oscillate to 
muon-neutrinos and contribute to the observable muon-signal. Calculations of limits on MSSM 
parameters from observed muon fluxes |Q, g, |9[ have to take this effect into account. In addition, 
present and future high-energy neutrino telescopes as AMANDA [jlO| , ANTARES |jn| , BAIKAL 



|12| and NESTOR [13| have large detecting volumes which make a search for cascades induced 



by electron- or tau-neutrinos promising. It should be mentioned that event rates due to tau- 



neutrinos have been previously calculated in |14] for a specific SUSY model 



In the next section we discuss some aspects of the MSSM phenomenology relevant for the 
production of tau-neutrinos. In section |3| we discuss the propagation of tau-neutrinos through 
the interior of the sun, taking into account scattering, absorption and regeneration as well as the 
possibility of neutrino mixing. We then show how Uj- <-5- v^ oscillations modify the expected muon 
flux from the sun and the earth. We examine the prospects of direct detection of tau-neutrinos. 
Following 1 15, 16 1 we choose model-independent, extreme cases of the branching ratios, as well 



as specific intermediate branching ratios motivated by the MSSM phenomenology. 

2 MSSM phenomenology 

The four neutralinos are a linear combination of the super-partners of the electroweak gauge 

bosons and the Higgs boson. In the basis {B,Wz,Kl,h%), the neutralino mass matrix is given 

by 

^ Ml —mzsm.9y/ cos (3 ttt.^ sin ^vy sin /? \ 

M2 rriz cos By/ cos (3 —rriz cos 9y/ sin /3 

mz sin 6v/ cos /? —mz cos 9w cos j3 — // 

y TTT,^ sin 0VF sin /? — tti^cos^vk sin/3 — /i J 

where Mi and M2 are the gaugino mass parameters, /i is the higgsino mass parameter and tan/3 
is the ratio of the Higgs vacuum expectation values. Here we use the convention and notation 



(2) 



of 1 17]. The eigenvectors of @) then represent the four neutralinos of which the lightest, XO) is 

Xo = N^iB + iVoaVFs + NqJi\ + N^Jil (3) 

We define the gaugino fraction to be R-^^ = Nq^ + Nq2 and refer to a neutralino as gaugino-like if 
R^ > 0.9, as higgsino-like if R^ < 0.1, and as mixed if the gaugino fraction is inbetween 0.1 and 
0.9. To obtain the distribution of branching ratios we have evaluated the cross-sections given 



in |17] for annihilation into cc, bb, tt,TT, W~^W~ and Z^Z^ in the approximation of WIMPs 
being at rest. Figure || shows, by means of different shading, the dominant branching ratios in 
the M2 — /U plane, with the trilinear mass chosen to be Ma_=1000 GeV, tan/3=3 and degenerate 
sfermion masses of 200 GeV. The GUT relation Mi = |M2 x tan'^ 9^ was assumed to hold. The 
mass contours of the lightest neutralino, xo, as well as the gaugino fraction (i?^=0.1, 0.9) are 
indicated by the full and dashed lines, respectively. 

A higgsino-like neutralino is already ruled out provided its relic density lies within the cos- 
mologically interesting region 0.1 < Qy^h"^ < 0.3 M. 
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Figure 1: Distribution of branching ratios for the final states cc^ bb,tt,TT,W~^W and Z^Z^. 
The final states with the largest branching ratio are indicated by different shading. The cross- 



sections |17| were evaluated in the limit v^ — > 0. We assumed Ma=WOO GeV, tan/3=3, degener- 



ate sfermion masses of 200 GeV and Mi = |M2 x tan Ow- Contours of equal M-^ are indicated 
by the full lines. Dashed lines distinguish regions of different gaugino/higgsino character. 



To estimate the order of magnitude of the different branching ratios, we evaluate the cross- 
sections for a few simple neutralino compositions. In particular we focus on the ratio B^^j^/Br- = 
a^jur-^ since, if the neutralino is gaugino-like and lighter than the top-quark, annihilation into 
hh and tt is responsible for most of the expected neutrino-signal. The final states bh/rT also 
produce extremely low/high fluxes of tau-neutrinos compared to the fluxes of muon-neutrinos. 

Generally, if annihilation happens via s-channel exchange of the Z^ and neutral Higgs bosons, 
A, and if mass threshold effects are neglected, the ratio is given by 



B^T 1 m 



2 



'TT 



B- ^ ctm'^' ^'^' 

^bb *-/ ^b 

where cj = 3 is the number of colors and m-r and m^ are the tau and b-quark masses, respectively. 

The mass dependency is a result of the helicity-suppressed annihilation into two fermions. Here, 

QCD correction to a^ becomes significant but can be absorbed into a "running mass" |17]. To 



give an example, for a neutralino of 80 GeV, mi) reduces to 3/4 of its tree-level value. We thus 
obtain a ratio of about Bj-y/B^j^ ~ 0.1 . 

If the neutralino would be a pure B-ino (A^'oi = 1 and A''o2 = Nq-^ = Nq4 = 0), the charge 
dependent t- and u-channel exchange of sfermions dominates the cross section, and we obtain 

B^ I m^ml / (T, - erf + ei y 

Bb-b Cfmlmf[{Ts-e,y + el) ' ^^ 

Here, T3 is the weak isospin (= —1/2) and e is the electrical charge (e,- = — 1 and Cb = —1/3). 
It was assumed that the slepton and squark masses, m^- and rriq, are horizontally degenerated. 
For mj = rriq and m-^ =80 GeV it follows that Btt/B^^ ~ 9. It should be noted that if the GUT 
relation Mi = |M2Xtan^ 9w holds and the neutralino is gaugino-like it usually also resembles the 
B-ino, but with a small higgsino contribution remaining. This higgsino contribution typically 
weakens relation (|5|). By scanning over a large range of MSSM parameters we find a more 
characteristic ratio for B-ino like neutralinos of Bj-y/B^ = 0{1). Note also, that this ratio 
increases for m^ < rriq as predicted in mSUGRA models ||18|| . 

If the neutralino is W-ino like (A'^02 = 1 and Nqi = Nq^ = Nq4 = 0) and for a neutralino mass 
range m]yz < ^^x ^ "^t' annihilation into gauge bosons generally dominates. For m^ < myy^z^ 
the annihilation into fermions is characterized by equation (^. It should be noted that even 
though the GUT relation Afi = tM2 x tan^ 9w prohibits the lightest neutralino to be W-ino like 



|17|, there are models with anomaly-mediated supersymmetry breaking which predict a W-ino 



like LSP |9| 

If the neutralino is of mixed composition, we can expect a large branching ratio for annihi- 
lation into W~^W~ and Z^Z^ [20|. Finally, if the neutralino is heavier than the top mass, the 



branching ratios for annihilation into a pair of top quarks is almost always favored. 

Summarizing, it can be said that the branching ratio for tt is highest when the lightest 
neutralino is primarily B-ino and not heavier than the top-quark. Thus it is this area of the 
MSSM parameter-space for which significantly higher fluxes of tau-neutrinos can be expected. 

3 Neutrino propagation: scattering, absorption, regeneration 
and oscillation 

Neutrinos in the energy range relevant for WIMP-annihilation propagate nearly unscattered 
through the earth. However, for the sun, its large size and high density in the core make it 



necessary to include neutrino reactions |14]. The main differences between tlie reactions of tau- 
and electron- or muon-neutrinos is that the charged current (CC) cross-section for tau-neutrinos 
is reduced due to kinematic reasons, and that a CC tau-neutrino reaction will regenerate a 
tau-neutrino in the prompt decay of the fully polarized tau. Thus the tau-neutrinos have to be 
treated differently. A particular interesting situation arises when neutrino mixing is considered. 



The general case of neutrino mixing in the presence of media was discussed in j21] and leads in 



case of forward-scattering to the well studied MS W-effect ||2^ . The propagation of oscillating 
neutrinos can be described by their density matrix |2^], pp-. 

Pp = -i[np, Pp] + / dp'{Wp>pPp> - Wppipp) (6) 

1 . r ^ 

-{pp,Ii}Ap+ / dpIrPp'IrRp'p 



'^flavor 

■ E 



where ilp is the momentum-dependent matrix of oscillation frequencies, Wpp' is the transition 
probability for NC-scattering of a neutrino of momentum p into a state p' , which is a product of 
the density and the composition-dependent cross-section. It is the projection operator onto the 
flavor i, At is the absorption rate of flavor i due to a CC-reaction and Rpip is the regeneration 
probability of a tau-neutrino. 

The first term in equation (|6|) is responsible for flavor-oscillations while the second term, 
since the cross-sections for neutral current (NC) reactions are to a good approximation flavor- 
independentj^, describes coherent scattering. It is similar to a gain-loss term in a usual Boltzmann- 
type collision integral. The third term accounts for absorption, which is flavor dependent due 
to the smaller CC cross-section of the tau-neutrino. However, the fourth and last term breaks 
the coherent development by adding pure tau-neutrinos due to the regeneration-mechanism de- 
scribed above. We thereby neglect the muon- or electron-neutrinos which are produced in 17 % 
of the tau-neutrino regeneration. 

Here, we will restrict ourself to two flavor Vr ^^ i/^ oscillations. For the energies under 
consideration, this can be justified as mixing of z^e •*-> i^^ is suppressed due to the low mass 
differences of '^ 10^^ eV^ derived from solar neutrino data |g] and due to the relatively short 
baseline of the high density core of the sun. Mixing of Ve ^^ t^t is also partially suppressed due 



to constraints on the effective mixing-angle, sin^ 29^^ < 0.1 |23]. 

Because of the strongly inhomogeneous density and composition of the sun, the analytic 
treatment of (y) becomes very involved and we thus treat the problem by means of a Monte 
Carlo calculation. We simulate the neutrino propagation by stochastic interactions along its 
way out of the sun. The cross-sections for NC- and CC-reactions were calculated using the MRS 



parton distribution function 1 24 ] . The density and composition of the sun as a function of the 
radius were taken from [p^] . The regeneration of the tau-neutrino by tau-decay was simulated 
according to [ 26 1 . The dynamics of equation (^ is incorporated in the following way: A coherent 



development of the oscillating neutrino-beam is assumed even when the neutrino undergoes a 
NC-reaction. In case of a CC-reaction the neutrino is absorbed and in the case of a tau-neutrino 
regenerated as a pure tau-neutrino with a lower energy. 

Figure y shows the average neutrino-energy after propagating from the core to the surface 
of the sun as a function of the injection energy. It can be seen that, due to tau-neutrino 



We therefore neglect the MSW-effect related case of Ve — e forward scattering. 
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Figure 2: Average neutrino-energy after propagating from the core to the surface of the sun, 
as a function of the injection energy. The injected neutrino beam is initiahy single flavored but 
after propagation represents the sum of the oscillating flavors. The full lines show the case of no 
oscillation while the dashed lines corresponds to full mixing and Am^ = 0.003 eV^. The upper 
two curves correspond to an initial tau-neutrino beam and the lower two represent an initial 
muon-neutrino beam. Electron-neutrino attenuation is similar to the case of muon-neutrino 
attenuation without oscillation. Absorbed neutrinos are included in the average with E^=0. 



regeneration, the tau-neutrino beam is considerably less attenuated than the muon- or electron- 
neutrino beam and that the difference between the mixing and no-mixing case is relatively 
small. The latter can be understood by comparing the typical absorption length, XabsorpUon-, 
to the oscillation length, XoscUiation, and the size of the scattering region dgcatter- Neutrino 
oscillations can have a strong influence if 

'^oscillation r^ '^absorption r^ "•scatter ■ \' ) 

For a neutrino-energy of 200 GeV, a total cross-section of 0.5 x 10^'^^ cm^ x £'i^(GeV), an 
average core density of 100 g/cm^ and a scattering region of dgcatter — 0.1 x Rq (which contains 
approximately 60% of the total column depth) we obtain 



^oscillation 



^absorbtion — -^ X Uscatter (o) 



Hence, it can be argued that above 200 GeV the influence of neutrino oscillations is reduced due 
to the long oscillation-length and below this energy, scattering becomes less important. This 
means that in a first approximation, propagation through the sun decouples from neutrino- 
oscillations and can therefore be treated separately. We note that a similar line of arguments 
applies to the influence of oscillation in the scenario of 4 neutrino flavors in which the fourth 
neutrino is a sterile neutrino. 



We conclude this section with a comment on the relevance of the MSW-effect |22], which 



has been discussed for the relevant energy range in |27|. Resonant flavor-conversion between the 
electron- and muon-neutrino is generally suppressed due to the high neutrino-energies consid- 
ered here. However, present constrains on the oscillation-parameter space allow resonant flavor 
conversion between the electron- and tau-neutrino if E^^ ^ < 75 GeV |27]. The conditions for the 
MSW-effect are met, if at all, in the outer layers of the sun, where the neutrinos practically not 
suffer from further attenuation. The above results are therefore unaffected, and we will neglect 
the MSW-effect in the further discussion. 

4 Event Rates 

To calculate event rates for a given detector, formula (|^) has to be folded with the detection 
probability of the neutrino. Here we will assume a detector of 1 km x 1 km x 1 km dimensions 
and, somewhat arbitrarily, set its energy threshold to Ey = 25 GeV. For a discussion of the effect 
of the energy threshold see H]. Since detailed muon rate calculations exist for non-oscillating 
muon-neutrinos [Q, ^ we will restrict ourself to the analysis of changes of the event rates due to 
oscillation or contained tau-events. 

We follow [^, |l^ and evaluate the extreme cases of branching ratios, Bp, for which signal 
expectations are either lowest or highest. If the neutralino is lighter than the top-quark, the 
hardest neutrino spectrum can be expected when the neutralino annihilates into pairs of taus, 
W^s or Z^s [28|. We thus assume an upper limit given by an exclusive annihilation into taus. 



A lower limit is given by annihilation into a pair of 6-quarks, with the branching ratio being 
almost always signiflcantly larger than for annihilation into c-quarks. If the neutralino is heavier 
than the top-quark the dominant annihilation channel will be the top quark or, if the neutralino 
is higgsino-like, it might also annihilate into pairs of W s or Z^s. We neglect final states which 
include gluons, since they produce generally a very soft neutrino spectrum which is not suitable 



for detection in large scale neutrino-detectors |29|. We will also omit final states containing 



Higgs bosons. The decay of the neutral Higgs is characterized by an equation similar to (^ or, 
if heavy enough, by annihilation into W s or Z^s. Annihilation into charged Higgs bosons can 



lead to an increase of the tau- neutrino flux, as was recently discussed in [3C]. 

The decay of heavy-quarks, gauge-bosons and taus was simulated with JETSET ||3^. Stop- 
ping of the b- and c-quark in the sun can become significant and we included it by sampling 



from a stopping spectrum developed in |14]. The distribution of WIMPs is assumed to be con- 
centrated at the core of the sun which is a valid approximation for the WIMP-masses considered 
here |jl3]. The neutrino propagation through the sun has been discussed above. The muon range 
in water is approximately calculated by 



L{E, 



'iJ^^ 



Pp 



a/ 13 



a/(3 



(9) 



with Q=0.0025 GeVg ^cm-^, /3 = 4.0 x 10 ^g ^cm^ and p being the density of the medium [| 
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Figure 3: Ratios of neutrino-induced muon fluxes with and without oscillations, (p^^osci/^'fi^ are 
shown for the sun (left) and the earth (right). The cases of extreme and intermediate branching 
ratios (assuming Bi^f, + Brr = 1) are shown by the full four lines, each one for a different Btt- 
Annihilation into W-bosons (top quarks) is represented by the dotted (dashed) line. Annihilation 
into Z-bosons leads to a similar ratio as annihilation into W-bosons. In the case of i/^ <-^ i^r 
oscillations, full mixing and Am^ = 0.003 eV^ has been assumed. The detector was assumed to 
have an energy threshold of Ei, = 25 GeV. 



^-)- v^ oscillations on the neutrino-induced muon flux for the 



Figure 13 shows the effect of v. 
sun and the earth in the extreme cases discussed above, as well as for some moderate cases, 
motivated in section |2[ The oscillation case is normalized to the no-oscillation case. When 
annihilating exclusively into a pair of taus, the muon flux increases due to oscillations up to 
a factor of 4-6 for the sun and 3-4 for the earth, which is explained by the higher and harder 



tau- neutrino flux |14|. This enhancement is only weakly dependent on the exact value of the 
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Figure 4: Ratios of CC contained tau-event rates, Rr, to muon-neutrino event rates (contained 
and semi-contained), R^, in the absence (upper part) and presence (lower part) of Ur ^^ u^ 
oscillations. Annihilation into fermions (Bt-t + -655 = 1) is represented by the four lines. The 
case of annihilation into W-bosons (top-quarks) is represented by the dotted (dashed) line. A 
detector of km x km x km dimensions with an energy threshold of £'1/ = 25 GeV has been 
assumed. 



branching ratio, and even a small branching ratio of B-^t = 0.1 leads to considerably higher 
muon fluxes. Only for a ratio of Br-/B^j^ < 0.05, a deficit due to neutrino oscillations can be 
expected. By scanning a large range of MSSM parameters as well as from equations (^) and (^), 
we find that such a small ratio is below of what is expected in most MSSM models. Above the 
top-mass threshold, annihilation into taus as final states can be neglected. Here, the difference 
between the two extremes B]yw=^ and Btt=l is small. Due to the long baseline in the case of 
the sun, oscillations lead to a small increase in the rate. For the earth, oscillations produce no 
visible effect, since for higher neutralino- masses the average oscillation length exceeds the radius 
of the earth considerably. Annihilation into Z^Z^ produces a rather similar ratio as annihilation 
into W+W-. 

Generally we find, that for a wide range of MSSM models, I'r *-^ i^^i oscillations lead to equal 
or larger muon fluxes. We note that this is a different conclusion than drawn in [^, where no 
increase but a generic decrease in rates due to neutrino-oscillations was found. 

The ratio of contained tau- neutrino event rates, Rr, to the muon- neutrino event rate, R^, is 
shown in figure ^ for the case of no oscillations (upper part) and the case of Uj- <-> z/^ oscillations 
(lower part). Here, only detection by CC-reactions has been considered. The muon event 
rate includes contained and semi-contained events, where the vertex is located outside of the 
detection volume. It was already shown in [Q] that the contribution of the semi-contained events 
is small for a detector of the large size assumed here. With increasing WIMP mass, the average 
muon range grows and therefore leads to a decrease in the ratio Rr/Rfj,- This is seen for the 
case of the earth, but for the sun, this effect is compensated by the larger attenuation of the 
muon-neutrinos within the sun. If no oscillation is present, the event rate due to tau-neutrinos 
is significantly larger than due to muon-neutrinos, as long as the branching ratio B-rr is large. 
In the case of Ur <-5- z/^ oscillations, higher tau- neutrino event rates are only preserved for WIMP 
annihilation in the earth. For the sun, the flux of tau- and muon neutrinos completely average 
out, and the difference in event rates are due to the lower cross-section of the tau neutrino and 
the additional semi-contained muon events. 

5 Discussion 

After reviewing briefly some of the MSSM phenomenology relevant for tau-neutrino production, 
we have discussed the propagation of WIMP-neutrinos through the sun and the earth under 
consideration of the presently most favored u-j- <-> v^ oscillation-parameters. Am? = 0.003 eV^ 
and sin 2^ = 1 [0]. A main result is that in a first approximation, scattering, absorption and 
regeneration of neutrinos in the sun can be effectively thought of as decoupled from oscillation- 
effects. Therefore the attenuation of the neutrino beam can be calculated first, and effects due 
to flavor-mixing can be applied later. 

We have shown that, Uj- ^^ v^ oscillations can alter the muon flux from WIMP-annihilation 
significantly. In particular if the lightest neutralino is B-ino like and lighter than the top- 
quark, we expect an increase of the muon fiux of up to a factor of 4-6. If annihilation happens 
primarily through gauge-bosons or the top-quark, oscillations will lead to only small changes in 
the muon-flux, as the disappearing of the original muon-neutrinos, f^ — > i/,-, is compensated by 
the appearance of original tau-neutrinos, z^t- — > i/^. A decrease due to i^r "^ t^fi oscillations is 
observed only if annihilation happens primarily into fermions and if Brr/Bbi, ^ 0.05, a ratio of 
branching ratios which is lower than expected for most of the MSSM parameter space. 

Because of the significantly higher tau-neutrino flux for parts of the MSSM parameter space, 

10 



contained tau-events are a promising signal if a large detection volume is available. However, 
the high tau-neutrino flux can be, relative to the muon-neutrino flux, considerably reduced if 
muon and tau-neutrinos mix. We omit electron-neutrinos from our discussion, since their flux 
at the origin is very similar to that of the muon-neutrinos. 

When comparing the detection potential for the different channels it is important to note that 
the background of atmospheric neutrinos sets limits to the WIMP-neutrino sensitivity, and that 
the backgrounds due to muon-, electron- or tau-neutrinos are significantly different in rate. In 
the relevant energy range, the flux of atmospheric muon-neutrinos is 10-20 times higher than the 



electron- neutrino flux |34| and unless f^ -^ v^ oscillations are present, the flux of tau-neutrinos 
is negligible [35|. Additionally there is also a small neutrino-background from the sun itself 
|36|. Experiment-related constraints might favor the detection of the CC muon-neutrino events 
which generally consist of long muon-tracks in contrast to tau- or electron-neutrinos producing 
cascade-like events. However, since contained cascade-like events deposit their energy within the 
detector we can expect a good energy resolution. This might add a new option to the rejection 
of atmospheric- neutrino events |37]. In case of a measurable WIMP signal rate, an energy 



spectrum of cascade-like events could be used to get a rough estimate of the WIMP-mass. By 
comparing the fluxes of muon- and tau-neutrinos it might become possible to obtain a hint on 
the composition of the WIMP. 
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